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Nonuniversality and finite dissipation in decaying magnetohydrodynamic turbulence 
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A model equation for the Reynolds number dependence of the dimensionless dissipation rate in 
freely decaying homogeneous magnetohydrodynamic turbulence in the absence of a mean magnetic 
field is derived from the real-space energy balance equation, leading to Ce ~ Ce,ao + CjR- + 

0{1/R?-)), where R- is a generalized Reynolds number. The constant Ce,oo describes the total 
energy transfer flux. This flux depends on magnetic and cross helicities, because these affect the 
nonlinear transfer of energy, suggesting that the value of Ce,ao is not universal. Direct numerical 
simulations were conducted on up to 2048® grid points, showing good agreement between data and 
the model. The model suggests that the magnitude of cosmological-scale magnetic fields is controlled 
by the values of the vector field correlations. The ideas introduced here can be used to derive similar 
model equations for other turbulent systems. 

PACS numbers: 47.65.-d, 52.30.Cv, 47.27.Jv, 47.27.Gs 


Magnetohydrodynamic (MHD) turbulence is present 
in many areas of physics, ranging from industrial appli¬ 
cations such as liquid metal technology to nuclear fusion 
and plasma physics, geo-, astro- and solar physics, and 
even cosmology. The numerous different MHD flow types 
that arise in different settings due to anisotropy, align¬ 
ment, different values of the diffusivities, to name only 
a few, lead to the question of universality in MHD tur¬ 
bulence, which has been the subject of intensive research 
by many groups [1-12]. The behavior of the (dimension¬ 
less) dissipation rate is connected to this problem, in the 
sense that correlation (alignment) of the different vector 
fields could influence the energy transfer across the scales 
[2, 13, 14], and thus possibly the amount of energy that 
is eventually dissipated at the small scales. 

For neutral fluids it has been known for a long time 
that the dimensionless dissipation rate in forced and 
freely decaying homogeneous isotropic turbulence tends 
to a constant with increasing Reynolds number. The 
first evidence for this was reported by Batchelor ]15] in 
1953, while the experimental results reviewed by Sreeni- 
vasan in 1984 [16], and subsequent experimental and nu¬ 
merical work by many groups, established the now well- 
known characteristic curve of the dimensionless dissipa¬ 
tion rate against Reynolds number: see [17-20] and ref¬ 
erences therein. For statistically steady isotropic turbu¬ 
lence, the theoretical explanation of this curve was re¬ 
cently found to be connected to the energy balance equa¬ 
tion for forced turbulent flows [19], where the asymptote 
describes the maximal inertial transfer flux in the limit 
of infinite Reynolds number. 

For freely decaying MHD, recent results suggest that 
the temporal maximum of the total dissipation tends to 
a constant value with increasing Reynolds number. The 
first evidence for this behavior in MHD was put forward 
in 2009 by Mininni and Pouquet [21] using results from 
direct numerical simulations (DNSs) of isotropic MHD 
turbulence. The temporal maximum of the total dissipa¬ 
tion rate e{t) became independent of Reynolds number 


at a Taylor-scale Reynolds number Rx (measured at the 
peak of e{t)) of about 200. 

Dallas and Alexakis [22] measured the dimensionless 
dissipation rate Cg from DNS data, where e was non- 
dimensionalized with respect to the initial values of the 
rms velocity U{t) and the integral length scale L{t) (here 
defined with respect to the total energy), for random ini¬ 
tial fields with strong correlations between the velocity 
field and the current density. The authors compared data 
with Ref. [21], and again it was found that Cg —)■ const. 
with increasing Reynolds number. Interestingly the ap¬ 
proach to the asymptote was slower than for the data of 
Ref. [21]. 

In this Letter we propose a model for the Reynolds 
number dependence of the dimensionless dissipation rate 
derived from the energy balance equation for MHD tur¬ 
bulence in terms of Elsasser fields [23], which predicts 
nonuniversal values of the dimensionless dissipation rate 
in the infinite Reynolds number limit. In order to com¬ 
pare the predictions of the model against data, we carried 
out a series of DNSs of decaying isotropic MHD turbu¬ 
lence. Firstly we explain the derivation of the model 
equation, then proceed to a description of our numerical 
simulations and subsequently compare the model to DNS 
results. We conclude with a discussion of the results and 
suggestions for further research. 

The equations describing incompressible decaying 
MHD flows are 

dtU = - VP — (tt • V)tt H— {\/X b) X b + nAu , (1) 

P P 

dtb = {b ■ V)m — (tt • V)b -1- ijAb , (2) 

V • M = 0 and V • b = 0 , (3) 

where u denotes the velocity field, b the magnetic induc¬ 
tion expressed in Alfven units, n the kinematic viscosity, 
r] the resistivity, P the pressure and p = 1 the density. 
For simplicity and in order to compare to results in the 
literature we consider the case of unit magnetic Prandtl 
number, that is Pm = vjr} = 1. 
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For freely decaying MHD turbulence the decay rate of 
the total energy Ed = —dtEtot equals the total dissipa¬ 
tion rate e, and the time evolution of the total energy 
is governed by the energy balance equation of MHD tur¬ 
bulence in real space, which is derived from the MHD 
equations (l)-(3). This suggests that the energy balance 
equation can be used in order to derive the Reynolds 
number dependence of the total dissipation rate. 

Since we are interested in the total dissipation e = 
^mag + Ekin, where ^mag = 2r] dk PEjnagik) , and 
£fcm = dk k'^Ekinik) {Emag{k) and Ekin{k) denot¬ 

ing magnetic and kinetic energy spectra), are the mag¬ 
netic and kinetic dissipation rates, respectively, we could 
take two approaches, either formulating the energy bal¬ 
ance in terms of the primary fields u and b or in terms 
of the Elsasser fields = u ± b. Since 

dt{\z^\^) =2dtEtot±2dtH, , (4) 

where El,, = {u ■ b) is the cross helicity, we can de¬ 
scribe the total dissipation either by the energy balance 
equations for {\z^\^) [23] or by the sum of the energy 
balance equations for Emag{t) = dk Emag{k) and 
EMn{t)= S^dkE^nik) [24, 25], 

This, however, is not the case if we are interested in 
the dimensionless dissipation rate. Unlike in hydrody¬ 
namics, there are several choices of scales with which to 
non-dimensionalize e(t) and thus with respect to which 
to define an MHD analogue to the Taylor surrogate ex¬ 
pression [15, 18]. For example U and L could be used, or 
the rms b field B and L or U and Lkin etc., or scales 
defined with respect to z^. The physical interpreta¬ 
tion is different for the different scaling quantities. Since 
the total dissipation must equal the total flux of energy 
passed through the scales by the kinetic and magnetic 
energy transfer terms, a scaling with U will be appropri¬ 
ate only for hydrodynamic transfer as this transfer term 
scales as U^/Lkin- All other transfer terms include b 
and u and thus should be scaled accordingly. This also 
precludes the most straightforward generalization of the 
Taylor surrogate, which would be a scaling of e with L 
and \/U'^ + B"^. A hydrodynamic transfer term would 
then be scaled partly with magnetic quantities, while the 
appropriate scaling should only involve kinetic quantities. 

Instead we propose to define the dimensionless dissipa¬ 
tion rate for MHD turbulence with respect to the Elsasser 
variables 


a = 


Ct + C: 



( 5 ) 


where L± = (Sn dk k Wz^\^))/{A dk (| 2 =^P)) 
are the integral scales defined with respect to z^, and 
denote the rms values of 2 :^ [26]. 

Using this definition we can now consistently non- 
dimensionalize the evolution equations of (| 2 ^P). For 


conciseness we outline the arguments for the (| 2 +p) case, 
since the {\z~\^) case proceeds analogously [27]. 

Following [23] the energy balance for (| 2 +p) reads for 
the case Pm = 1 

+ ('''UrBtt) . ( 6 ) 

where (r) and B^^ (r) are the longitudinal third- 

order correlation function and the second-order longi¬ 
tudinal structure function of the Elsasser fields, respec¬ 
tively. The definitions of the these functions can be found 
in the Supplemental Material [27]. Using (4) one can ex¬ 
press the LHS of (6) in terms of e(t) and dtHc- 

If we now introduce the nondimensional variable a = 
r/L+ [3] and nondimensionalize equation (6) with re¬ 
spect to 2 ^ and L+ as proposed in the definition of 
in eq. (5), we obtain 


r<+ — " 


22+^2- 


-u2 _ 




-u2 _ 


-dtH, 


■ r] idc 


L+z~ 2(7'^ 


71 ) + + 

^4 a ^LL 

^ ^(7 . n 


( 7 ) 


In this way we arrive at a consistent scaling for each 
transfer term in (6) with the appropriate quantity, as the 
function scales with 2 +^ 2 “. 

Since the inverse of the coefficient in front of the dissi¬ 
pative term is similar to a Reynolds number, we introduce 
the generalized large-scale Reynolds number 


R- = 


22 Lj^ 
n + ri 


( 8 ) 


hence (7) suggests a dependence of C+ on IjR-. How¬ 
ever, the structure and correlation functions and the cross 
helicity flux also depend on Reynolds number. 

For conciseness we introduce dimensionless versions of 
all terms present on the RHS of (7), such that 

= z+^z-g+-+{a,t), (9) 

Bttir,t) = z+\++ia,t), (10) 

dtBt+ir,t) = ^-^^^E+{a,t), ( 11 ) 

dMt) = (^^fcG+(t), (12) 


which leads to a dimensionless version of the (| 2 +p) en¬ 
ergy balance equation for freely decaying MHD turbu¬ 
lence 


Ct = 




3cr4 +_ + \ 
) 


R- Cr4 




^-F+ + G+ 
4 


( 13 ) 
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After non-dimensionalization the highest derivative in 
the differential equation is multiplied with the small pa¬ 
rameter 1 /-R_, suggesting that this can be viewed as a 
singular perturbation problem [28]; and thus we consider 
asymptotic expansions of the dimensionless functions in 
inverse powers of i?_ [19, 29]. 

The formal asymptotic series of a generic function / 
(used for conciseness in place of the functions on the RHS 
of (7)) up to second order in l/i?_ reads 

f = fo + + (^-^2 + 0{RJ') . (14) 

After substitution of the expansions into (13) and fol¬ 
lowing the same steps for the evolution equation for 
{\z~\^), we arrive at model equations for Cf and Cj 

cf = + ^ + + 0{R-^) , (15) 


temporal dependences of e.g. should be separable in 
the inertial range, that is ~ (£“''(t)r)“ for some 

value a, and dtB^^ ~ Q;£+(t)““^9t£+r“. At the peak of 
dissipation = 0 , 

and we obtain Ff^{tpgak) = 0. As the terms GJ which 
describe the flux of cross helicity in the infinite Reynolds 
number limit, cancel the corresponding contribution from 
the transfer terms [27], the asymptotes C^ao describe the 
flux of total energy provided the model (15) is applied at 

ipeak• 

Due to selective decay, that is the faster decay of the 
total energy compared to He and H^ag [14], one could 
perhaps expect dtHe to be small compared to £ in the 
inhnite Reynolds number limit in most situations. In 
this case we obtain GJ ~ 0 and 

Cte^itpeak) = ~ ’ ( 21 ) 


up to third order in l/i?=p, where we defined the coeffi¬ 
cients and 


Q'ioo = -§ f ± ^0^ - (16) 


G± = 




±T± 


= 




,4 I o ,±± _ 5l 


±± 




a dahfi — 


- ± G± , 


(17) 
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(18) 


in order to write (13) in a more concise way. Using = 
{L-/L+){z^/z~)R- to define 


G = 


1 

2 



A_ z+ 
L+ z~ 


c- 



(19) 


{D is defined analoguously), finally one obtains for the 
dimensionless dissipation rate Gg 

Gg = Gg^oo + —I" + 0{R_^) . (20) 


Since the time dependence of the various quantities 
in this problem has been suppressed for conciseness, we 
stress that (20) is time dependent, including the Reynolds 
number i?_. A normalization using initial values of z^ 
and would have resulted in a dependence of Gg(t) on 
initial values of i?_, which only describe the initial condi¬ 
tions and not the evolved flow for which Gg is measured. 

At the peak of e{t) the additional terms should 
in fact vanish for constant flux of cross helicity (that is, 
d^Hc = 0), since in the infinite Reynolds number limit 
the second-order structure function will have its inertial 
range form at all scales. By self-similarity the spatial and 


which recovers the inertial-range scaling results of 
Ref. [23] and reduces to Kolmogorov’s 4/5th law for 

b = 0. 

Since Gg^oo is a measure of the flux of total energy 
across different scales in the inertial range, differences 
for the value of this asymptote should be expected for 
systems with different initial values for the ideal invari¬ 
ants He and magnetic helicity Hmag = (a • h), where a is 
the vector potential b — W x a. In case of H^nag 7 ^ Oj the 
value of Gg^oo should be less than for Hmag = 0 due to 
a more pronounced reverse energy transfer in the helical 
case [13] [30], the result of which is less forward transfer 
and thus a smaller value of the flux of total energy. For 
He 7 ^ 0 we expect Gg_oo to be smaller than for He = 0 , 
since alignment of u and b weakens the coupling of the 
two fields in the induction equation, which leads to less 
transfer of magnetic energy across different scales and 
presumably also less transfer of kinetic to magnetic en¬ 
ergy. In short, one should expect nonuniversal values of 

Cs,oo- 

Before we compare the model equation with DNS data 
and address this question of nonuniversality numerically, 
we briefly outline our numerical method. Equations (1)- 
(3) are solved numerically in a periodic box of length 
Lbox = 27r using a fully de-aliased pseudospectral MHD 
code [31, 32]. All simulations satisfy kmaxVmag,km ^ 
where rimag,kin are the magnetic and kinetic Kolmogorov 
scales, respectively. We do not impose a background 
magnetic field, and both the initial magnetic and ve¬ 
locity fields are random Gaussian with zero mean, with 
initial magnetic and kinetic energy spectra of the form 
Emag,kin{k) ~ exp(-fc2/(2A:o)^), where /cq > 5 and 

further simulation details are specified in Table 1 of 
[27]. The initial relative magnetic helicity is Pmag{k) = 
kHmag{k)/2Emag{k) = 1 for all runs of series H and zero 
for the runs labelled NH. The initial relative cross he¬ 
licity was yOc(O) = i7c(0)/(lM(0)||6(0)|) = 0 for runs of 
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the H and NH series and Pc(0) = 0.6 for series CH06H 
and CH06NH, while initial magnetic and kinetic ener¬ 
gies were in equipartition. All spectral quantities have 
been shell- and ensemble-averaged, with ensemble sizes 
restricted by computational resources to up to 10 runs 
per ensemble. The total dissipation rate e was measured 
at its maximum. 

Figure 1 shows hts of the model equation to DNS data 
for datasets that differ in the initial value of Hmag and 
He- As can be seen, the model fits the data very well. 
For the series H runs and for i?_ > 70 it is sufficient to 
consider terms of first order in i?_, while for the series 
NH the first-order approximation is valid for R_ > 100. 
The cross-helical CH06H runs gave consistently lower 
values of compared to the series H runs, while lit¬ 
tle difference was observed between series CH06NH and 
NH. The asymptotes were Ce,ao = 0.241 ± 0.008 for 
the H series, = 0.265 ± 0.013 for the NH se¬ 
ries, = 0.193 ± 0.006 for the CH06H series and 

^,00 = 0.268 ± 0.005 for the CH06NH series. 

As predicted by the qualitative theoretical arguments 
outlined before, the measurements show that the asymp¬ 
tote calculated from the nonhelical runs is larger than for 
the helical case, as can be seen in Fig. 1. The asymptotes 
of the series H and NH do not lie within one standard er¬ 
ror of one another. Simulations carried out with ^ 0 
suggest little difference in Cg for magnetic fields with 
initially zero magnetic helicity. For initially helical mag¬ 
netic fields Ce is further quenched li He ^ 0. In view 
of nonuniversality, an even larger variance of Ce^oo can 
be expected once other parameters such as external forc¬ 
ing, plasma /3, Pm, etc., are taken into account. Here 
we have restricted ourselves to nonuniversality caused by 
different values of vector field correlations. 

In summary, a definition for the dimensionless dissi¬ 
pation rate for MHD turbulence has been proposed, 
where e was non-dimensionalized with respect to the 
Elsasser fields instead of the rms velocity. For this def¬ 
inition of Cg and the case of unit Prandtl number we 
derived a model for the dependence of Cg on a gener¬ 
alized Reynolds number i?_. The model predicts that 
Ce —t const with increasing i?_, in analogy to hydro¬ 
dynamics, and the asymptote is a measure of the total 
energy transfer flux. The model was compared to DNS 
data for datasets which differ in their initial values of 
magnetic and cross helicities. At moderate to high i?_, 
we found good agreement to data with the model only 
using terms up to first order in 1 /i?_. However, at low 
terms of second order in i?_ cannot be neglected, in 
fact these terms improve the fit specifically at low i?_. 
This is expected from adding another term in the expan¬ 
sion and thus provides further justiheation of the validity 
of eq. (20). 

As predicted, the values of the respective asymptotes 
from the datasets differ, suggesting a dependence of 
on different values of the helicities, and thus a connection 



FIG. 1. (Color online) The solid and dotted and dash-dotted 
lines show (20) fitted to helical, non-helical and cross-helical 
DNS data, respectively. The red (grey) lines refer to fits using 
the first-order model equation, the black lines use the model 
equation up to second order in 1/R-. As can be seen, the 
respective asymptotes differ for the data sets. 


to the question of universality in MHD turbulence. This 
presents an interesting point for further research concern¬ 
ing the influence of other vector held correlations on the 
dissipation rate. Other questions concern the general¬ 
ization of this approach to more general MHD systems 
such as hows with magnetic Prandtl numbers Pm ^ 1, 
compressive huctuations, and to the presence of a back¬ 
ground magnetic held, as well as to turbulent systems 
where the how carries other quantities such as tempera¬ 
ture or pollutants; and also the application to decaying 
hydrodynamic turbulence [20]. In the most general case 
in plasmas there will be a mean magnetic held, which 
leads to spectral anisotropy and the breakdown of the 
conservation of magnetic helicity [33] and thus might in¬ 
troduce several difficulties to be overcome when general¬ 
izing this method, as the spectral hux will then depend 
on the direction of the mean held [3, 34] and a more 
generalized description and role for the magnetic helicity 
would be needed. 

Our model shows that different degrees of correlation 
in a turbulent plasma control the amount of energy that 
can effectively be transferred into the smallest scales. It 
could have several possible applications, e. g. for heating 
rates in the solar wind, especially as high values of the 
cross helicity inhibit such transfer to some extent. For 
situations where one is interested in sustaining a mag¬ 
netic held over long times, thus trying to minimize dis¬ 
sipative effects, one could estimate from (16)-(17) what 
type of correlations produce not only a low asymptotic 
value of the dissipation rate but also a fast approach to 
this asymptote. This would have relevance to cosmo¬ 
logical and astrophysical [35] magnetic helds as well as 






5 


terrestial plasmas, such as in a tokamak reactor. Our 
results suggest that in cosmology, where a topical prob¬ 
lem is the origin of large-scale magnetic fields, it is not 
only a nonzero value of magnetic helicity, but perhaps 
also the parameter range of other correlations such as 
the cross and kinetic helicities, that facilitate the pres¬ 
ence of long-time magnetic fields. Moreover, this raises 
questions about the possible generation mechanisms for 
cosmological magnetic fields leading to different correla¬ 
tions between the vector fields such that they can sustain 
long evolution times. 
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In this Supplemental Material we include definitions of the structure and correlation functions 
used in the Letter. For simulations with zero cross helicity we give numerical verification that the 
two model equations for C'f' and Cf in fact lead to identical results, as expected. Moreover, we 
provide a table of simulation details with information on the resolution of our simulations, Reynolds 
numbers and measured values of the dimensionless dissipation rate. 


DEFINITIONS OF THE STRUCTURE AND CORRELATION FUNCTIONS 

In order to keep this material self-consistent, we include here the definitions of the Elsasser structure and correlation 
functions used in the Letter. The third-order correlation and structure functions l(^) 

and the second-order structure functions are defined as follows: 



= {{5LZ^{r)f5LZ^{r)), 

( 1 ) 



( 2 ) 



( 3 ) 

Bti{r) 

= {{5LZ^{r)f), 

( 4 ) 


where vl = v ■ rjr denotes the longitudinal component of a vector field v, that is its component parallel to the 
displacement vector r, and 




[v(a; -h r) 



( 5 ) 


its longitudinal increment. 

The transfer term in the energy balance equations for 2 :^ are written in the Letter in terms of the function J (r) 
for reasons of conciseness. It is common in the literature to use the relation 




( 6 ) 


in order to express the transfer terms through the functions ^ (r) and J (r) instead [1]. 


VALIDATION FOR DIFFERENT SCALINGS FOR ZERO CROSS HELICITY 


The definition of Cg as in eq. (5) of the Letter was proposed in order to arrive at a consistent scaling for the transfer 
terms in the energy balance equations for , since the functions which describe the transfer terms scale 

as This led to separate model equations for (7+ and C~. In the Letter, we outlined the derivation of the 

model equation for C+ from rescaling the evolution equation for (| 2 ;+p) leading to 


cy = c,:„ + |L + ^ + o(Ry). 


D+ 


( 7 ) 


In the same way a model equation for 


C7 = 



( 8 ) 


is derived through nondimensionalizing the energy balance equation for z . The flux terms in this equation are given 
in terms of the function C'£'^^(r) (or equivalently and B~f^~f^{r)), and the dissipative term in terms of 



2 




FIG. 1. The solid and dotted lines show the model equations fitted to helical and non-helical DNS data, respectively. Figure 
(a) shows and Fig. (b) . In both figures, the red lines refer to fits using the first-order model equations, the black lines 

use the model equations up to second order in 1/R^±. By comparing Figs, (a) and (b) it can clearly be seen that both scalings 
give the same results. 




FIG. 2. The solid lines show the model equations for C^, and fitted to data for (a) series GH06H and (b) series CH06NH 
with initial relative cross helicity Pc(0) = 0.6. The red crosses refer to additional runs with 0.2 ^ Pc(0) 0.8. 


For zero cross helicity one should expect C+ = , since 

(|z+|2) = 2Etot + 2Hc = 2Etot and {\z-\^) = 2Etot - 2H, = 2Etot ■ (9) 

Therefore all quantities defined with respect to the rms fields z'^ and z~ should be the same, and thus for zero cross 
helicity it is possible to use, Cf or C~ representatively for Cg. 

Figure 1 shows separate fits of the model equations for (7+ and C~ (obtained by the two nondimensionalizations) 
to DNS data for zero cross helicity. The measured asymptotes are identical and the curves are practically indistin¬ 
guishable from each other and from those shown in Fig. (1) of the Letter, which was obtained by fitting eq. (20) of 
the Letter to the data. This provides numerical verification for the representative use of (7+ (or C~) to describe Cg. 

However, for nonzero cross helicity the difference in the rms values and precludes the use of (7+ (or C~) 
to represent the dimensionless dissipation rate C^, and both (7+ and C~ have to be taken into account separately, 
which is reflected in the definition of the dimensionless dissipation rate in eq. (5) of the Letter. Figure 2 shows data 
for nonzero cross helicity He, that is runs of series CH06H and CH06NH, as well as the additional series CHH and 
CHNH. Interestingly, very little difference between (7+ and C~ is observed for series CH06H. For series CH06NH, 
note that the asymptote Ce,oo is the same, no matter if extrapolated using or i?+. 
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SIMULATION SPECIFICATIONS 


Run id 


^maxT^mag 

R- 

Rl 

R\ 

V 

ko 

# 

c. 

a 

pc(0) 

HI 

128^ 

1.30 

33.37 

25.28 

14.87 

0.009 

5 

10 

0.756 

0.008 

0 

H2 

256^ 

2.42 

37.77 

27.81 

15.85 

0.008 

5 

10 

0.704 

0.007 

0 

H3 

512^ 

1.38 

50.81 

35.08 

18.55 

0.002 

15 

10 

0.608 

0.001 

0 

H4 

256® 

1.80 

61.14 

40.63 

20.34 

0.005 

5 

10 

0.569 

0.006 

0 

H5 

256® 

1.59 

76.72 

48.73 

23.11 

0.004 

5 

10 

0.510 

0.005 

0 

H6 

1024® 

1.38 

89.32 

55.51 

25.76 

0.00075 

23 

10 

0.4589 

0.0003 

0 

H7 

256® 

1.29 

102.53 

60.65 

26.91 

0.003 

5 

10 

0.450 

0.004 

0 

H8 

512® 

2.33 

123.17 

69.40 

29.67 

0.0025 

5 

10 

0.419 

0.003 

0 

H9 

512® 

2.01 

154.67 

83.06 

33.84 

0.002 

5 

10 

0.384 

0.003 

0 

HIO 

512® 

1.45 

255.89 

123.97 

45.21 

0.0012 

5 

10 

0.320 

0.004 

0 

HU 

528® 

1.31 

308.69 

143.71 

50.18 

0.001 

5 

10 

0.310 

0.004 

0 

H12 

1024® 

2.03 

441.25 

194.38 

61.39 

0.0007 

5 

5 

0.281 

0.002 

0 

H13 

1032® 

1.38 

771.34 

309.08 

82.97 

0.0004 

5 

5 

0.268 

0.001 

0 

H14 

1024® 

1.24 

885.05 

358.72 

88.76 

0.00035 

5 

5 

0.265 

0.002 

0 

H15 

2048® 

1.35 

2042.52 

724.71 

136.25 

0.00015 

5 

1 

0.250 

- 

0 

CH06H1 

512® 

2.17 

124.89 

108.81 

49.88 

0.002 

5 

1 

0.380 

- 

0.6 

CH06H2 

512® 

1.57 

207.61 

171.87 

68.57 

0.0012 

5 

5 

0.309 

0.002 

0.6 

CH06H3 

1024® 

2.21 

351.52 

277.21 

95.31 

0.0007 

5 

1 

0.260 

- 

0.6 

CH06H4 

1024® 

1.76 

491.50 

380.70 

116.85 

0.0005 

5 

1 

0.236 

- 

0.6 

CH06H5 

1024® 

1.37 

696.19 

523.08 

132.48 

0.00035 

5 

1 

0.231 

- 

0.6 

CHHl 

512® 

1.46 

254.64 

129.83 

47.55 

0.0012 

5 

5 

0.315 

0.002 

0.2 

CHH2 

512® 

1.50 

240.35 

147.68 

55.25 

0.0012 

5 

5 

0.311 

0.003 

0.4 

CHH3 

512® 

1.74 

149.28 

195.27 

90.60 

0.0012 

5 

5 

0.323 

0.03 

0.8 

NHl 

256® 

1.51 

55.57 

53.89 

25.57 

0.004 

5 

10 

0.587 

0.005 

0 

NH2 

256® 

1.26 

71.51 

68.60 

30.11 

0.003 

5 

10 

0.530 

0.004 

0 

NH3 

512® 

1.86 

103.41 

96.69 

37.68 

0.002 

5 

10 

0.468 

0.004 

0 

NH4 

512® 

1.51 

133.14 

122.51 

43.94 

0.0015 

5 

10 

0.431 

0.004 

0 

NH5 

512® 

1.29 

161.35 

151.76 

50.73 

0.0012 

5 

10 

0.394 

0.004 

0 

NH6 

1024® 

2.28 

192.40 

168.28 

54.44 

0.001 

5 

5 

0.358 

0.002 

0 

NH7 

1024® 

1.76 

259.58 

232.10 

65.42 

0.0007 

5 

5 

0.358 

0.002 

0 

NH8 

1024® 

1.40 

354.30 

301.71 

76.73 

0.0005 

5 

5 

0.323 

0.002 

0 

NH9 

2048® 

1.15 

1071.44 

823.58 

134.73 

0.00015 

5 

1 

0.279 

- 

0 

CH06NH1 

512® 

2.02 

94.39 

113.02 

49.29 

0.002 

5 

1 

0.482 

- 

0.6 

CH06NH2 

512® 

1.41 

148.86 

174.61 

65.48 

0.0012 

5 

5 

0.417 

0.003 

0.6 

CH06NH3 

1024® 

1.93 

242.06 

272.85 

87.50 

0.0007 

5 

1 

0.365 

- 

0.6 

CH06NH4 

1024® 

1.52 

325.62 

365.54 

104.45 

0.0005 

5 

1 

0.341 

- 

0.6 

CH06NH5 

1024® 

1.16 

450.01 

515.23 

127.09 

0.00035 

5 

1 

0.313 

- 

0.6 

CHNHl 

512® 

1.31 

168.31 

152.84 

51.40 

0.0012 

5 

5 

0.401 

0.003 

0.2 

CHNH2 

512® 

1.34 

162.27 

157.48 

55.52 

0.0012 

5 

5 

0.402 

0.006 

0.4 

CHNH3 

512® 

1.58 

116.71 

204.78 

87.52 

0.0012 

5 

5 

0.456 

0.002 

0.8 


TABLE I. Specifications of simulations. H refers to an initially helical magnetic field, NH to an initially non-helical magnetic 
field. Rl denotes the integral-scale Reynolds number, R\ the Taylor-scale Reynolds number, R- the generalized Reynolds 
number as in given in eq. (8) of the Letter, rj the magnetic resistivity, ko the peak wavenumber of the initial energy spectra, 
kmax the largest resolved wavenumber, r/mag the Kolmogrov microscale associated with the magnetic field at the peak ot total 
dissipation, # the ensemble size, the dimensionless total dissipation rate defined in eq. (5) of the Letter, a the standard 
error on Ce and pc(0) the initial relative cross helicity. All Reynolds numbers are measured at the peak of total dissipation. 
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ENERGY AND CROSS-HELICITY FLUXES 


In the Letter it is pointed out that the asymptotes Cf' describe the total energy flux, as the contribution of the 
cross-helicity flux to the Elsasser flux is cancelled by the respective terms . Here we provide some further details. 
The asymptotes were given in the Letter as 


e.oo ^4 


dr, (3a^ ±T±\ _ 5 


-9o 




which reduced at the peak of dissipation to 




The dimensional version of this equation is 




( 10 ) 


(11) 


( 12 ) 


where we assume that the function has its inertial range form corresponding to ■ The function 

(13) 


can be written in terms of and 


= 4 (^LL,i - 2C'i0j = - (((SLZ^(r)ySLZ^(r)} - 2{zf{x)zf{x)z];{x + r))) , 
which can be written in terms of the primary helds u and b as 

C'lUL = \ - 2{z^{x)z^{x)z]^{x + r))) 

= 7x((<5lu(0)^ - QbL{xfuL{x + r)) T HnSLiHr))^ - QuL{xfbL{x + r)) , 


43 


43 


(14) 


see e.g. [1], and the two terms on the last line of (14) are related to the flux terms in the evolution equations of the 
total energy and the cross helicity [1]. 

Now going back to (12), we can write in the inertial range 


dr ( 3r^ 


e = —- 




a ^-4 


— [{{6Lu{r)f - &bL{xfuL{x -b r)) =F {{5L{b{r)f - QuL{xfbL{x + r))] ) ± dtH^ 


= et ± ec ± dtHc = St , 


(15) 


where St is the flux of total energy and sc the cross-helicity flux, which must equal —dtHc in the inertial range for 
freely decaying MHD turbulence, thus we see that the contribution from the third-order correlator ^ resulting in 
Sc is cancelled by dtHc, or, after nondimensionalization, the cross helicity flux scL±/[{z^)‘^z~] is cancelled by GJ. 
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